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We investigate the pp —* pK + A reaction within an effective Lagrangian model where the contribu- 
tions to the amplitudes are taken into account within the tree level. The initial interaction between 
the two nucleons is modeled by the exchange of 7T, p, u> and a mesons and the AK + production 
proceeds via the excitation of the iV*(1650), iV*(1710), iV*(1720) baryonic resonances. The pa- 
rameters of the model at the nucleon-nucleon-meson vertices are determined by fitting the elastic 
nucleon-nucleon scattering with an effective interaction based on the exchange of these four mesons, 
while those at the resonance vertices are calculated from the known decay widths of the resonances 
as well as the vector meson dominance model. Available experimental data is described well by 
this approach. The one-pion-exchange diagram dominates the production process at both higher 
and lower beam energies. The p and uj meson exchanges make negligible contributions. However, 
the cr-exchange processes contribute substantially to the total cross sections at lower beam energies. 
The excitation of the A* (1710) and A* (1650) resonances dominate this reaction at beam momenta 
above and below 3 GeV/c respectively. The interference among the amplitudes of various resonance 
excitation processes is significant. For beam energies very close to the K + production threshold the 
hyperon-proton final state interaction effects are quite important. The data is selective about the 
model used to describe the low energy scattering of the two final state baryons. 
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I. INTRODUCTION 

In recent years there has been a considerable amount 
of interest in the study of the associated production re- 
action p + p — ► p + K + + A. This is expected to provide 
information on the manifestation of quantum chromody- 
namics (QCD) in the non-perturbative regime of energies 
larger than those of the low energy pion physics where 
the low energy theorem and partial conservation of axial 
current (PCAC) constraints provide a useful insight into 
the relevant physics Q|. The strangeness quantum num- 
ber introduced by this reaction leads to new degrees of 
freedom into this domain which are expected to probe 
the admixture of ss quark pairs in the nuclcon wave 
function || and also the hyperon-nucleon and hyperon- 
strange meson interactions . 

The elementary nucleon-nucleon-strange meson pro- 
duction cross sections are the most important ingredients 
in the transport model studies of the if + -meson produc- 
tion in the nucleus- nucleus collisions, which provide infor- 
mation on not only the initial collision dynamics but also 
the nuclear equation of state at high density . Fur- 

thermore, the enhancement in the strangeness produc- 
tion has been proposed as a signature for the formation 
of the quark-gluon plasma in high energy nucleus- nucleus 
collisions US 

The experimental data on the pp — > pK + A reaction is 
very scarce. The measurements performed in late 1960s 
and 1970s provide total cross sections for this reaction at 
beam momenta larger than 2.80 GeV/c (see e.g. [p^|). 
With the advent of the high-duty proton-synchrotron 
COSY at the Forschungszentrum, Jiilich, it has become 
possible to perform systematic studies of the associated 



strangeness production at beam momenta very close to 
the reaction threshold (2.340 GeV/c). The first round 
of experiments at COSY have already added |l6| 12 new 
data points to the data base. At near threshold beam 
energies the final state interaction effects among the out- 
going particles are significant. Therefore, the new set of 
data are expected to probe also the hyperon-nucleon and 
hyperon-strange meson interactions along with the mech- 
anism of the strangeness production in proton-proton col- 
lisions. 

The existing theoretical studies of this reaction are 
based either on a single boson (tt or K meson) exchange 
mechanism |l7]-|2Cj or a resonance model p2|-p4|. In the 
first method, the K + production is assumed to take place 
essentially through the exchange of one intermediate pion 
or if -meson; the excitation of any intermediate nucleon 
resonance is not considered. The AT-meson exchange am- 
plitudes were found to dominate [fl7|j2(J the production 
cross sections. However, the relative sign of the pion and 
AT-meson exchange amplitudes was not fixed in this ap- 
proach pC|| . Furthermore, it has been argued that the 
existing high energy data can be well reproduced consid- 
ering only the single pion-exchange process Jll||l9| since 
the contribution of the A'-meson exchange amplitude can 
be compensated by various parameters of the model. 

In the resonance model [^]J of the strangeness produc- 
tion in pp collisions, the 7r, ry, and p-meson exchanges are 
included and the K + meson production proceeds via the 
excitation of the N* (1650) , N* (1710) and N* (1720) reso- 
nances [^2|,0. However, the terms in the total amplitude 
involving the interference of various resonance contribu- 
tions are neglected in these calculations. Moreover, the 
parameters of the NNir and NNp vertices were taken 
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from the Bonn nucleon-nucleon potential which may not 
be adequate at higher beam energies as these have been 
determined by fitting the NN scattering data below the 
NNir production threshold. At the same time, the finite 
life-time of the p meson has not been taken into account 
while calculating the relevant coupling constants from the 
experimental branching ratios. 

In this paper, we investigate the associated K + pro- 
duction in the proton-proton collisions in the framework 
of an effective Lagrangian approach |25|-|29[|, following 
and extending our previous study [ |30| , |3l[ on 7r° and 7r + 
production. Initial interaction between two incoming nu- 
cleons is modeled by an effective Lagrangian which is 
based on the exchange of the 7r, p, uj and a mesons. The 
coupling constants at the nucleon-nucleon-meson vertices 
are determined by directly fitting the T-matrices of the 
nucleon-nucleon (NN) scattering in the relevant energy 
region Q . The effective-Lagrangian uses the pseudovec- 
tor (PV) coupling for the nuclcon-nucleon-pion vertex 
(unlike the resonance model P2[), and thus incorpo- 
rates the low energy theorems [|34|] of current algebra 
and the hypothesis of partially conserved axial-vector 
current (PC AC). The K + production proceeds via ex- 
citation of the AT*(1650), AT* (1710) and A*(1720) inter- 
mediate baryonic resonance states which have apprecia- 
ble branching ratios for the decay into the K + A channel. 
The interference terms between various intermediate res- 
onance states are included which marks a major differ- 
ence between our work and the resonance model p2| . To 
describe the recent near threshold data, the final state 
interaction between the outgoing particles is included 
within the framework of the Watson-Migdal theory (3^] . 

The remainder of this paper is organized in the fol- 
lowing way. Section II contains details of our theoretical 
approach. Section III comprises the results of our analy- 
sis and their critical discussion. The summary and con- 
clusions of our work are presented in section IV. Finally 
some technical details are given in appendix A. 



A. Nucleon-nucleon-meson vertex 

The parameters for JVJV vertices are determined by fit- 
ting the NN elastic scattering T matrix with an effective 
NN interaction based on the 7T, p, u and a meson ex- 
changes. The effective meson- JVJV Lagrangians are 
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FIG. 1. Feynman diagrams for K + A production in pp 
collisions. The diagram on the left shows the direct process 
while that on the right the exchange one. 



C-NNp 

£-NNlo 
£nN<j 



9NNtt 



fan 



XT ■ p^^ N . 

-gNNu^N 



2itin 



(2.1) 

(2.2) 
(2.3) 
(2.4) 



In Eqs. (2) - (3) cr M „ is defined as 



II. FORMALISM 

We consider the tree- level structure (Fig. 1) of the am- 
plitudes for the associated K + A production in proton- 
proton collisions, which proceeds via the excitation of 
the JV*(1650)(±~), A*(1710)(± + ) and A*(1720)(f + ) in- 
termediate resonances. To evaluate these amplitudes 
within the effective Lagrangian approach, one needs to 
know the effective Lagrangians (and the coupling con- 
stants appearing therein) at (i) nucleon-nucleon-meson, 
(ii) resonance-nucleon-meson, and (iii) resonance- K + - 
hyperon vertices. These are discussed in the following 
subsections. 
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We have used the notations and conventions of Bjorken 
and Drell |||]. In Eq. (1) to at denotes the nucleon mass. 
It should be noted that we have used a PV coupling 
for the NNtt vertex. Since we use these Lagrangians 
to directly model the T-matrix, we have also included 
a nuclcon-nucleon-axial-vector-isovector vertex, with the 
effective Lagrangian given by 
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where A represents the axial- vector meson field. This 
term is introduced because in the limit of large axial 
meson mass (titia) it cures the unphysical behavior in 
the angular distribution of JVJV scattering caused by the 
contact term in the one-pion exchange amplitude |33[ , if 
9nna is chosen to be 
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with very large mjv) mA- fn appearing in Eq. (7) is 
related to g NNn as f n = {^^)m n . 

It should be noted that the contact term of the coordi- 
nate space potential corresponding to one pion exchange 
term, is effectively switched off by the repulsive hard core 
of the nucleon-nucleon interaction. However, in the effec- 
tive Lagrangian description this term has to be explicitly 
subtracted in order to avoid the unphysical behavior of 
the NN elastic scattering cross sections. This is achieved 
by the inclusion of a term corresponding to the exchange 
of a axial- vector-isovector meson as described above. 

We introduce, at each interaction vertex, the form fac- 
tor 
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where qi and m, are the four momentum and mass of 
the ith exchanged meson respectively. The form factors 
suppress the contributions of high momenta and the pa- 
rameter Ai, which governs the range of suppression, can 
be related to the hadron size. Since NN elastic scatter- 
ing cross sections decrease gradually with the beam en- 
ergy (beyond certain value), we take energy dependent 
meson-nucleon coupling constants of the following form 



g(Vs) = g exp(-£y/s), 



(2.9) 



in order to reproduce these data in the entire range 
of beam energies. The parameters, go, A and i were 
determined by fitting to the elastic proton-proton and 
proton-neutron scattering data at the beam energies in 
the range of 400 MeV to 4.0 GeV |||§f|. It may be 
noted that this procedure also fixes the sign of the ef- 
fective Lagrangians (Eq. (l)-(4) and (6)). The val- 
ues of various parameters are shown in Table 1, which 
are the same as those used in the calculations of the 
pion production in proton-proton collisions [|30| , j3l"| . Thus 
we ensure that the NN elastic scattering channel re- 
mains the same in the description of various inelas- 
tic channels within this approach, as it should be. 

TABLE I. Coupling constants for the AW-meson vertices 
used in the calculations 



Meson 


<?74tt 


£ 


A 

(GeV) 


mass 
(GeV) 




12.562 


0.1133 


1.005 


0.138 


a 


2.340 


0.1070 


1.952 


0.550 




46.035 


0.0985 


0.984 


0.783 


P 


0.317 


0.1800 


1.607 


0.770 


k p = 6.033, k^ = 0.0 











B. Resonance- nucleon- meson vertex 

As the A-hyperon has zero isospin, only isospin-1/2 
nucleon resonances are allowed. Below 2 GeV center of 
mass (cm.) energy, only three resonances, TV* (1650), 
iV*(1710), and 2V*(172Q), have significant decay branch- 
ing ratios (3-11%, 5-25% and 1-15% respectively |D) 
into the K + A channel. In this work only these three res- 
onances have been considered. The JV*(1700) resonance 
having very small (and uncertain) branching ratio for the 
decay to this channel has been excluded. 

Since all the three resonances can couple to the meson- 
nucleon channel considered in the previous section, we re- 
quire the effective Lagrangians for all the four resonance- 
nucleon-meson vertices corresponding to all the included 
resonances. For the coupling of the spin- 1/2- resonances 
to pion we again have the choice of PS or PV couplings. 
The corresponding effective Lagrangians can be written 
as MM 
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where M — (mjv ± m^r), with the upper sign for even 
parity and lower sign for odd parity resonance. The op- 
erators T, T M , are given by, 
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for resonances of even and odd parities, respectively. We 
have performed calculations with both of these couplings. 
The effective Lagrangians for the coupling of resonances 
to other mesons are similar to those given by Eq. (2)- (4), 
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The operators V and T ^ v are, 
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for resonances of even and odd parities, respectively. 

The even parity isospin-1/2 iV*(1720) resonance is a 
spin-3/2 nucleon excited state. We have used the fol- 
lowing effective Lagrangians for vertices involving this 
resonance p9|,|36|] 
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Here, ^ M is the TV* (1720) vector spinor. It should be 
remarked here that an operator Q ai _ l (z) = g au — |(1 + 
2z)7 Q 7 AI has also been included in the vector spinor ver- 
tex in Refs. |2^|j{|. This operator describes the off-shell 
admixture of the spin- 1/2 fields f}7|| . The choice of the 
off-shell parameter z is arbitrary and it is treated as a 
free parameter to be determined by fitting to the data. 
This operator can be easily introduced in Eqs. (19)-(22) 
which will introduce four additional free parameters in 
our calculations. We however, work with the Lagrangians 
as given in Eqs. (19)-(22), which are identical to those 
given in 1 29 3^] for z = 0.5. 



C. Resonance- hyperon-strange meson vertex 

For vertices involving spin- 1/2 resonances, there is 
again the PS and PV coupling option. In principle one 
can select a linear combination of both and fit the PS/PV 
ratio to the data. However, to minimize the number of 
parameters we choose either PS or PV coupling at a time. 
The effective Lagrangians for the N*AK + vertex is writ- 
ten in the following way Q , 

For spin- 1/2 resonance, 
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ratios for the decay of the resonances to the correspond- 
ing channels. 

The partial width for the decay of a resonance (in its 
rest frame) of mass M/v* into a meson of mass m m and 
energy E m and a nucleon is written in terms of the Lorenz 
invariant matrix element M as 



dT 



(27T) 



2M N , 

m,N d 3 pN 
X E N (2ir) 3 



\M\ 2 5\P N * - Pm -p N ) 



d 3 p„ 



(2ir) 3 2E n 



(2.26) 



In case of the meson (in the decay channel) having a finite 
life time for the decay to another channel (e.g p — > tttt ), 
an integration over the phase-space for this decay must 
be included 138 



(i). N*Nir vertex 

For the spin- 1/2 resonance the N*Nir decay width, 
with the PS coupling, is given by 
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where 
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For spin-3/2 resonance the N*Nir decay width is 
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where M' = m^* ± m\, with the upper sign for even 
parity and lower sign for odd parity resonance. 
For spin-3/2 resonance, 
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D. Coupling constants for resonances 
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The plus and minus sign in Eqs (27) corresponds to odd 
and even parity resonances respectively, while in Eqs. 
(28) and (32) reverse is the case. 

(ii). N*Np vertex 

The partial decay width of each resonance for the decay 
into nucleon and two pions via the p meson is given by 



The resonance couplings are determined from the ex- T(tojv* 
perimentally observed quantities such as the branching 
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In this equation the spectral function S(m) is defined as, 
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The value of is taken to be 150 MeV. The p -> irir 

decay four- moment a p p7r7r are 

[to 2 — 4to 2 ] [m 2 ] 



Am 2 
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In Eq. (33) T*(m) is defined in the following way, 
For spin- 1/2 even parity resonance 



r* (to) = — 

V / A — 



i_ r 9N' 1/2 N P 

4ir \ 2m,N 

4(E* N + E m )(p™) 2 + 3(E* N - to n )to 2 



TUN' 



<J(p c ™) 2 +rn 



2 

N' 



E* N - 



E rn = Vfc m ) 2 +TO 2 , 



(2.38) 
(2.39) 
(2.40) 



where p^™ is given in the same way as Eq. (29) with to w 
replaced by m. 

For spin- 1/2 odd parity resonance 
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For spin-3/2 even parity resonance 
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(iii). N*Nlu vertex 



Since the resonances considered in this study have no 
known branching ratios for the decay into the Nlo chan- 
nel, we determine the coupling constants for the N*Nw 



vertices by the strict vector meson dominance (VMD) 
hypothesis , which is based essentially on the assump- 
tion that the coupling of photons on hadrons takes place 
through a vector meson. 

The N*N"f partial widths are given as following: 
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In these equations, qf = {{m 2 v , — to 2 v )/2totv*)- The value 
of /ijv* is determined by fitting to the N"f partial width 
of each resonance which is given in terms of the helicity 
amplitudes A 1 / 2 and A 3 / 2 by |55) 
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where J is the resonance spin, /in* is written as the ratio 
of the couplings at N*lu and W7 vertices as 
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Using above equations together with the experimental 
helicity amplitudes the values of the coupling constants 
for the N*Nuo vertices can be determined. We have used 
<7^7 = 17 in our calculations. 

(iv). N*Na vertex 

As the sigma meson is, most of the time, a resonance of 
two pions |42}| in the S-state, the coupling constants for 
the N*Na vertices are determined from the branching 
ratios of the decay of the resonances into N(ttttY =0 . We, 
however, reduce the experimental values of these ratios 
by two third to account for the fact that this correlated 
state provides only about 2/3 of the total 2vr-exchange. 
The expressions for the partial widths are similar to those 
given by Eqs. (27)-(32). 

(iv). N*AK+ vertex 

The coupling constants for the N*AK + vertex are de- 
termined from the experimental branching ratios for the 
N* — > A K + decay. The expressions for the decay widths 
are similar to those given by Eqs. (27) - (32). 

We assume that the off-shell dependence of the NN* 
vertices are determined solely by multiplying the ver- 
tex constants by the form factors, which have the dipole 
form 131131 
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The resonance properties used in the calculations of 
the decay widths are given in Table 2, where the resulting 
coupling constants and the adopted values of the cut-off 
parameters (Xf N ) are also given. It may be noted that 
we have fixed the latter to one value in order to minimize 
the number of free parameters. 

It should however, be stressed that the branching ra- 
tios determine only the square of the corresponding cou- 
pling constants; thus their signs remain uncertain in 
this method. Predictions from independent calculations 
(e.<7 the quark model) can, however, be used to con- 
strain these signs. The magnitude as well signs of the 
coupling constants for the N*Nir, N*AK, N*Np, and 
N* N(irir) s - wave vertices were determined by Feuster and 
Mosel J5^] and Manley and Saleski [Q in their analysis of 
the pion-nucleon data involving the final states ttN, ttttN, 
rjN, and KA. Predictions for some of these quantities are 
also given in the constituent quark model calculations of 
Capstick and Roberts (4^| . Guided by the results of these 
studies, we have chosen the positive sign for the coupling 
constants for these vertices. Unfortunately, the quark 
model calculations for the N*Nw vertices are still sparse 
and an unambiguous prediction for the signs of the corre- 
sponding coupling constants may not be possible at this 
stage [J46|. Nevertheless, we have chosen a positive sign 
for the coupling constants for these vertices as well. 

TABLE II. Coupling constants and cut-off parameters for 
the N* iV-meson and iV*-hyperon-meson vertices used in the 
calculations 



E. Propagators 

We require the propagators for various mesons and nu- 
cleon resonances in the calculation of the amplitudes. 
The propagators for pion, p meson and axial-vector 
mesons are given by 
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In Eq. (51), the mass of the axial meson is taken to be 
very large (188 GeV), as the corresponding amplitude is 
that of the contact term. The propagators for to and a 
mesons are similar to those given by Eqs. (50) and (49) 
respectively. 

The propagators for spin-1/2 and spin-3/2 resonances 
are 
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Resonance 


Width 

(GeV) 


Decay 
channel 


Adopted 
value of the 
branching 
ratio 




cut-off 
(GeV) 


AT* (1710) 


0.100 


Ntt 


0.150 


0.0863 


0.850 






Np 


0.150 


1.3653 


0.850 






Ncu 




0.1189 


0.850 






No 


0.170 


0.0361 


0.850 






AK 


0.150 


2.9761 


0.850 


AT* (1720) 


0.150 


Ntt 


0.100 


0.0023 


0.850 






Np 


0.700 


90.637 


0.850 






Nlu 




22.810 


0.850 






Ng 


0.120 


0.1926 


0.850 






AK 


0.080 


0.0817 


0.850 


AT* (1650) 


0.150 


Nn 


0.700 


0.0521 


0.850 






Np 


0.08 


0.5447 


0.850 






Nlu 




0.2582 


0.850 






No 


0.025 


0.2882 


0.850 






AK 


0.070 


0.0485 


0.850 



In Eqs. (52) - (53), Ym* is the total width of the res- 
onance which is introduced in the denominator term 
(p 2 — mfy.) to account for the fact that the resonances 
are not the stable particles; they have a finite life time 
for the decay into various channels. Tjv* is a function 
of the center of mass momentum of the decay channel, 
and it is taken to be the sum of the widths for pion and 
rho decay (the other decay channels are considered only 
implicitly by adding their branching ratios to that of the 
pion channel) 
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Tn*->n p is calculated according to Eq. (33). Pat^a^ is 
taken to be 
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where £ is the orbital angular momentum of the reso- 
nance, is as defined in Eq. (29) and is given 
by the same equation with uin* replaced by p of Eqs. 
(52) - (53). To is taken to be the total on-shell width 
of the resonance minus the corresponding width for the 
nucleon-p meson decay channel. 



G 



F. Amplitudes and cross sections 

After having established the effective Lagrangians, 
coupling constants and form of the propagators, we can 
now proceed to calculate the amplitudes for various dia- 
grams associated with the pp — > pAK + reaction. These 
amplitudes can be written by following the well known 
Feynman rules |fi"7f and calculated numerically by follow- 
ing e.g. the techniques discussed in |3(J. The isospin part 
is treated separately. This gives rise to a constant fac- 
tor for each graph, which is unity for the reaction under 
study. It should be noted that the signs of various ampli- 
tudes are fixed by those of the effective Lagrangian den- 
sities, coupling constants and propagators as described 
above. These signs are not allowed to change anywhere 
in the calculations. 

The general formula for the invariant cross section of 
the p + p — > p + A + K + reaction is written as 



da = 



1 



2^[{ Pl -p 2 y-m%] (2tt) ! 



-5\P f - P^Afi 



n 

a=l 



d 3 p a 
E a ' 



(2.56) 



where Afi represents the total amplitude, P, and Pf the 
sum of all the momenta in the initial and final states, 
respectively, and p a the momenta of the three particles 
in the final state. The corresponding cross sections in the 
laboratory or center of mass systems can be written from 
this equation by imposing the relevant conditions. 



G. Final state interaction 

For describing the data for the pp — > pAK + , reaction at 
beam energies very close to the threshold, consideration 
of the final state interaction (FSI) among the three out 
going particles is important. As there exists no theory of 
the FSI effects in the presence of three strongly interact- 
ing particles, we follow here an approximate scheme in 
line exactly with Watson-Migdal theory of FSI In 
this approach the energy dependence of the cross section 
due to FSI is separated from that of the primary produc- 
tion amplitude. This method has been applied earlier to 
study the low momentum behavior of the rj meson fig] ] 
and pion spectra [ pO| , pl|5T measured in proton- proton 
collisions. We write for the total amplitude 



A fi = M f i(pp -> pAK+) ■ T ff , 



(2.57) 



where Mfi(pp — > pAK + ) is the primary production am- 
plitude as discussed above, while Tff describes the re- 
scattering among the final particles which goes to unity 
in the limit of no FSI. The latter is taken to be the coher- 
ent sum of the two-body on-mass-shell elastic scattering 
amplitudes of the particles involved in the final channel. 



Tff 



(2.58) 



where U represents the two-body on-shell elastic scat- 
tering amplitude (of the interacting particles pair j — k) 
in the three-body space with the ith particle being the 
spectator. £i and qi denote the partial wave and relative 
momentum of the j — k particle pair respectively. 

An assumption inherent in the approximation given by 
Eq. (65) is that the reaction takes place over a small re- 
gion of space, a condition fulfilled rather well in near 
threshold reactions involving heavy mesons. This al- 
lows us to express the amplitudes h in terms of the in- 
verse of the Jost function, Je^qi) ^2 31 1. In the analy- 
sis presented in this paper we assume £i = 0, for all 
the three pairs j — k and use the modified Cini-Fubini- 



Stanghellini formula 1 53 1 for the effective range expansion 
of the phase-shift (5oi) of the relevant pair 

Cl qiCotS 0i + 2<?;7 ? h{rj) = (1/a,) + (l/2)r ^ 2 , (2.59) 

to calculate the corresponding Jost function. It may be 
noted that in case of the pair j — k involving uncharged 
particle(s), the second terms on the left hand side of Eq. 
(67) vanishes and Cq goes to unity. In this equation 
roi and a,i are the effective range and scattering length 
parameters respectively for the j — k interacting pair, r] 
is the corresponding Coulomb parameter and 



2 2tt?7 

e 27rr) _ I > 



^-J n(n 2 + 7] 2 ) 

71—1 K 1 



0.5772- ln{rj). 



In this case we have 



= (UK))- 1 = 



(if + ^V J2 



V< + (r c j2)qf - i Qi 
where a is given by 

a=(l/r^)[l + (l + 2r c Jad 1/2 }, 



and af and are defined as 



2qiV h(r))) 



1 Oi 



r 2 ' 



(2.60) 
(2.61) 

(2.62) 

(2.63) 

(2.64) 
(2.65) 



It may be noted that for large qi, the amplitude ti goes 
to unity which is to be expected. The extrapolation of 
the scattering amplitude for the off-shell effects can be 
achieved by means of a monopole form factor [ pp[ . For a 
detailed discussion of the off-shell effects we refer to |54j . 

The factorization of the total amplitude into those of 
the FSI and primary production (Eq. (65)), enables 
one to pursue the diagrammatic approach for the lat- 
ter within an effective Lagrangian model and investigate 
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the role of various meson exchanges and resonances in 
describing the reaction. Moreover, in this way the FSI 
among all the three outgoing particles can be included. 
Although the meson-baryon interactions are weak, they 
can still be influential through interference. 

The parameters a and ro are very poorly known for 
the if + -nucleon and i^ + -hyperon systems as the corre- 
sponding scattering data are scarce at low energies. On 
the other hand, for the hyperon-nucleon system several 
sets of values for these parameters have been given by the 
Bonn-Julich [ pq ] and Nijmegen |5(| groups from their re- 
spective A — p interaction models. There is quite some 
variation in the values given in these sets. For the K + -p 
and K + -A systems we have adopted the values given in 
a recent effective Lagrangian model analysis of the avail- 
able K A and Kp data by Feuster and Mosel |37| . In any 
case, the cross sections are insensitive to the FSI effects 
in these channels. On the other hand, these effects are 
very important in A — p channel, and we have performed 
calculation of the corresponding FSI effects with all the 
sets of these parameters given by Bonn-Julich and Ni- 
jmegen groups (given in Table 3) in order to see if the 
results are sensitive to various models. More details will 
be given in the next section. 



We notice that the measured cross sections are repro- 
duced reasonably well by our calculations (solid line) for 
all the beam energies except for two lowest points. The 
FSI effects which are not included in these calculations, 
reduce the discrepancy between the experimental data 
and calculations at these beam momenta. This point is 
further discussed in the next subsection. 

TABLE III. Scattering length (a) and effective range (ro) 
parameters for the AiV scattering derived from models A, A, 
B and B of the Jiilich-Bonn group jiij and D, F and NSC 
of the Nijmegen group 



model 


a (singlet) 
(fm) 


ro (singlet) 
(fm) 


a(triplet) 
(fm) 


ro (triplet) 
(fm) 


A 


1.56 


1.43 


1.59 


3.16 


A 


2.04 


0.64 


1.33 


3.91 


B 


0.56 


7.77 


1.91 


2.43 


B 


0.40 


12.28 


2.12 


2.57 


D 


1.90 


3.72 


1.96 


3.24 


F 


2.29 


3.17 


1.88 


3.36 


NSC 


2.78 


2.88 


1.41 


3.11 



III. RESULTS AND DISCUSSION 

The theoretical approach presented in the previous sec- 
tion has been used to study the available data on the p + 
p — > p + A + K + reaction for beam energies ranging from 
just above the production threshold to about 10 GeV. 
In the results shown below, we have used PS couplings 
for both N*Nir and N*AK + vertices involving spin-1/2 
resonances of even and odd parities. However, calcu- 
lations have also been performed with the corresponding 
PV couplings. The cross sections calculated with this op- 
tion for the resonance-hyperon-kaon vertex deviate very 
little from those obtained with the corresponding PS cou- 
plings. However, the PV coupling for the N*Nir vertex 
leads to noticeably different results as is discussed below. 

A. Cross section data for beam energy above 2 GeV 

In Fig. 2 we show the comparison of our calculations 
with the experimental data for the total cross section for 
this reaction as a function of beam momentum for inci- 
dent energies above 2 GeV. In this figure we have investi- 
gated the role of various meson exchange processes in de- 
scribing the total cross section. The dotted, dashed, long- 
dashed and dashed-dotted curves represent the contribu- 
tions of 7r, p, to and a meson exchanges respectively. The 
contribution of the heavy axial meson exchange is not 
shown in this figure as it is negligibly small. The coherent 
sum of all the meson-exchange processes is shown by the 
solid line. The experimental points are taken from p5|. 



We note that the pion exchange graphs dominate the 
production process for all the energies. In comparison 
to this, the contributions of p and ui meson exchanges 
are almost insignificant. The p-meson exchange, which 
is a convenient way of taking into account the P wave 
part of the correlated two-pion exchange (CTPE) pro- 
cess, is rather week even in the low energy NN scatter- 
ing p2[ . With increasing projectile energy its contribu- 
tion decreases further. On the other hand, the a me- 
son exchange, which models the CTPE in the 7T7T S-wave 
and provides about 2/3 of this exchange in the low en- 
ergy NN interaction, plays a relatively more important 
role. This observation has also been made in the case of 
NN -> NNtt reaction g^gQ. Thus, the a meson 
exchange provides an efficient means of mediating the 
large momentum mismatch involved in the meson pro- 
duction reactions in NN collisions, particularly at lower 
beam momenta. 

The relative importance of the contributions of each 
intermediate resonance to the pp —> pAK + reaction is 
studied in Fig. 3, where the contributions of iV*(1650), 
V*(1710) and V*(1720) resonances to the energy depen- 
dence of the total cross section are shown by dotted, 
dashed and dashed-dotted lines respectively. Their co- 
herent sum is depicted by the solid line. It is clear that 
the contributions from the iV*(1710) and V*(1650) res- 
onances dominate the total cross section at beam mo- 
menta above and below 3 GeV/c respectively. Moreover, 
the interference terms of the amplitudes corresponding 
to various resonances are quite important. This result is 
in sharp contrast to the resonance model calculations of 
Refs. ElLp2Lp4|, where these terms were ignored. It must 
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again be emphasized that we have absolutely no freedom 
in choosing the relative signs of the interference terms. 
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FIG. 2. The total cross section for thep + p — » p + K + + A 
reaction as a function of the beam momentum. The dotted, 
dashed, long-dashed and dashed-dotted curves represent the 
contributions of n, p, u) and a meson exchanges respectively. 
Their coherent sum is shown by the solid line. The experi- 
mental data are from ]l5| . 

Looking in Table 2, one might naively expect the dom- 
inance of the TV* (1710) resonance everywhere as the cou- 
pling constants for the N*AK + and N*Nn vertices for 
the A''* (1710) resonance are about an order of magni- 
tude larger than those for TV* (1720) and iV*(1650) res- 
onances. In fact, the relative importance of various res- 
onances is determined by the dynamics of the reaction 
mechanism. As the beam energy rises above the K + pro- 
duction threshold, the excitation of the resonance lowest 
in energy is more probable in the beginning. However, 
with increasing beam energy the excitation of the higher 
energy resonances starts playing more and more impor- 
tant role. 

As mentioned earlier, the use of the PV coupling for 
the N*AK vertices (involving spin-1/2 even and odd par- 
ity resonances) makes insignificant changes in the cross 
sections. However, there is a clear preference for the PS 
coupling at the N*Nir vertices. This is shown in Fig. 
4, where the ratio of the total cross section obtained by 
using the PV (<Jpv) and PS (aps) couplings for these 
vertices, is shown as a function of the beam momentum. 
It is seen that apv is larger than aps at higher beam mo- 
menta while at lower ones the reverse is true. Clearly PS 
coupling for the N*Nir vertex provides a better descrip- 
tion of the beam energy dependence of the total cross 
section for the pp — > pK + A reaction. 
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FIG. 3. Contributions of N*(1650) (dotted line), N*(1710) 
(dashed line) and N*(1720) (dashed-dotted line) baryonic res- 
onances to the total cross section for the p + p — > p + K + + A 
reaction as a function of beam momentum. Their coherent 
sum is shown by the solid line. 



B. Cross section data for beam energies below 2 
GeV 

In Fig. 5 we compare the results of our calculations 
(with FSI effects included) with the recent data jl6| for 
the pp — > pK + A reactions at beam energies very close to 
the kaon production threshold. In this figure the total 
cross section is shown as a function of the excess energy 
(e) = y/s - ttln - rriK+ - toa, where y/s is the invariant 
mass. The FSI effects were included by following the 
procedure outlined in section 2.G. We have chosen |57| a 
= 0.065 + i0.040, r = -15.930 - i8.252 and a = -0.214, 
ro = —0.331 for the K + A and K + p systems respectively 
in all the calculations shown below. For the A — p system 
all the seven sets of the parameters as shown in Table 3 
were used. 

In Fig. 5a, the results obtained with the parameters 
sets of models A (dotted line), A (solid line), B (dashed 
line) and B (long-dashed-dotted line) of the Bonn-Julich 
group |55| are shown. It can be noted that all the four 
models provide similar results for the total cross sections 
at e ~ 150 MeV. However, at lower values of e, cross 
section calculated with models A, and A are larger than 
those of model B and B by a factor of about 2-3. More- 
over, there is a difference of a factor of more than 2 be- 
tween the results obtained with model A and A itself 
with the latter providing a better overall agreement with 
the data. We also show in this figure the results obtained 
without including the FSI effects (dashed-dotted line). It 
can be noted that the FSI effects are quite important in 
order to describe the experimental data. 
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FIG. 4. Ratio of the total cross section calculated with 
pseudovector and pseudoscalar couplings for the N*Nir ver- 
tex corresponding to spin-1/2 (even and odd parity) resonance 
for the same reaction as in Fig. 2, as a function of beam mo- 
mentum. 

The results obtained with models D, E and NSC of 
the Nijmegen group [[36| are shown in Fig. 5b. These 
three models produce almost identical results for all the 
values of e. However, while the data at the higher ex- 
cess energies are reproduced by all the three model quite 
well, they under-predict the cross sections at lower e by 
a factor of about 3. Therefore, while all the models of 
A — p interaction (considered in this work) provide an 
equally good description of the total cross section data 
at higher values of the excess energy, a difference of fac- 
tors of 2-3 occurs between their predictions at lower val- 
ues of e. Thus, the near threshold AK + production data 
in proton-proton collisions are sensitive to the S-wave A- 
nucleon interaction and may be used to distinguish be- 
tween various models proposed in the literature to de- 
scribe this interaction. We note that model A of the 
Bonn-Jiilich group provides the best overall description 
of the data, which has been used to account for the A — p 
FSI effects in all the calculations discussed subsequently. 

The individual contributions of various nucleon res- 
onances to the total cross section of the pp — > pAK + 
reaction near the production threshold is shown in Fig. 
6, as a function of the excess energy. In contrast to the 
situation at higher beam energies (pi a b > 3 GeV/c), the 
cross section is dominated by the A^*(1650) resonance 
excitation. This is in line with the observations made in 
Rcfs. [§f,§ij]. It may, however, be noted that in Ref. [§4j, 
the FSI effects have not been included and no comparison 
with the data at near threshold energies is shown. 
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FIG. 5. The total cross section for the p + p — > p + K + + A 
reaction very close to the K + production threshold as a func- 
tion of the excess energy (defined in the text). The FSI effects 
are included by using the scattering length (a) and effective 
range (ro) parameters for the K + — A and K + — p systems 
taken from the Ref. j57j] and those for the A — p system from 
the sets given by Julich-Bonn |35| and Nijmegen |56j groups. 
In the upper part (a) results obtained with the A — p pa- 
rameters of models A (dotted line), A (solid line), B (dashed 
line) and B (long-dashed line) of the former group are shown, 
while in the lower part (b) those of models D, F and NSC 
of the latter group are depicted. Results of the three mod- 
els of the Nijmegen group are indistinguishable from each 
other. In the upper part (a) results with no FSI effects 
are shown by dashed-dotted line. The experimental data are 
taken from 

Since AT* (1650) is the lowest energy baryonic reso- 
nance having an appreciable branching ratio for the de- 
cay to the AK + channel, its dominance in this reaction 
at beam energies near the kaon production threshold is to 
be expected. The contributions of other two resonances 
(AT* (1710) and A^*(1720)) are several orders of magni- 
tude smaller, therefore, the resonance-resonance interfer- 
ence terms are also very small. Thus, near threshold 
energies, this reaction proceeds preferentially via excita- 
tion of the AT* (1650) resonance. It may be noted that in 
Ref. [p4j , the FSI effects have not been included in these 
calculations. 

In Fig. 7, we show the contributions of various meson 
exchanges to this reaction at near threshold beam ener- 
gies. Various curves have the same meaning as in Fig. 
2. The one pion exchange graphs dominate the reaction 
in this energy regime as well. On the other hand, the 
individual contributions of the p and ui meson exchange 
processes are negligible. However, those of the a meson 
exchange are substantial in this energy regime. Thus, 
like near threshold pion production in proton-proton col- 
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lisions, the heavy scalar meson exchange plays an impor- 
tant role in this case too. It should however, be noted 
that the interference terms of various meson exchange 
processes are not negligible; contributions of various ex- 
change processes simply do not add up to the total cross 
section obtained by the coherent addition of various am- 
plitudes. 
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FIG. 6. Contributions of N*(1650) (dotted line), N*(1710) 
(dashed line) and N*(1720) (dashed-dotted line) baryonic res- 
onances to the total cross section for the same reaction as in 
Fig. 5, as a function of the excess energy. Their coherent sum 
is shown by the solid line. The FSI effects are included with 
a and ro parameters of the K + — p and K + — A systems being 
the same as those in Fig. 5 and those for the A — p system 
being taken from model A of the Jiilich-Bonn group [^5| . 



IV. SUMMARY AND CONCLUSIONS 

We investigated the associated K + A production in the 
proton-proton collisions at energies ranging from near 
threshold to about 10 GeV. This reaction is of inter- 
est as it provides the prospect of testing QCD in the 
non-perturbative domain at energies larger than the pion 
mass. In this paper our goal has been to investigate this 
reaction within an effective Lagrangian model which is 
proven to be very successful in describing the pion pro- 
duction in NN collisions. Most of the parameters of this 
model are fixed by fitting to the elastic NN T-matrix; 
this restricts the freedom of varying the parameters of 
the model to provide a fit to the data. 

The reaction proceeds via the excitation of the 
JV*(1650), iV*(1710) and A*(1720) intermediate nucleon 
resonant states. The coupling constants at vertices in- 
volving resonances have been determined from the ex- 
perimental branching ratios of their decay into various 
relevant channels. Unlike the NNt: vertex, there is no 
compelling reason to choose the pseudovector (PV) form 



for the N*AK + and N*Nn couplings (involving spin-1/2 
resonances of even and odd parities) and we investigated 
both the PV and pseudoscalar (PS) couplings at these 
vertices. To describe the data at the near threshold beam 
energies, we have included the FSI effects among the out- 
going particles by following Watson-Migdal theory which 
has been used before successfully to describe the NNij 
and NNtt reactions in the similar energy regimes. 
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FIG. 7. Contributions of n (dotted line), p (dashed line), to 
(long-dashed line) and a (dashed-dotted line) meson exchange 
processes to the total cross section for the same reaction as 
shown in Fig. 6, as a function of the excess energy. Their 
coherent sum is shown by the solid line. The FSI effects are 
included in the same way as in Fig. 6. 

With the same set of parameters, the model is able to 
provide a good description of the data at higher as well 
as near threshold beam energies. The one-pion-exchangc 
processes make the dominant contribution to the cross 
section in the entire energy regime. The individual con- 
tributions of the p and uj meson exchange diagrams is 
very small every where. Although, the interference terms 
of their amplitudes with those of other meson exchanges 
may still be noticeable. On the other hand, the a ex- 
change makes a relatively larger contribution at lower 
beam energies, confirming the earlier observation that 
the heavy scalar meson provides a means of mediating the 
large momentum transfer in near threshold TV TV-meson 
production processes. 

While at beam momenta larger than 3 GeV/c, the 
reaction proceeds predominantly via excitation of the 
iV*(1710) resonance, the process gets maximum contri- 
bution from the A* (1650) resonance at lower beam ener- 
gies. A very striking feature of our results is that the in- 
terference among various resonance contributions is very 
significant. Therefore, in the calculations of this reaction 
these terms should not be ignored. The near threshold 
data clearly favors the excitation of the iV*(1650) reso- 
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nance. Therefore, this reaction, in this energy regime, 
provides an ideal means of investigating the properties of 
this baryonic resonance. 

Unlike the NNir vertex where there is a clear prefer- 
ence for the PV coupling, as seen in the NNir data, the 
present reaction does not distinguish between PS and PV 
couplings at the N*AK + vertex involving spin-1/2 even 
or odd parity resonance. However, the PS coupling at 
the N*Nir vertex is clearly favored by the data. 

The near threshold data may be selective about the 
model describing the low energy A-nucleon scattering. 
Calculations of the FSI effects performed with the scat- 
tering length and effective range parameters of the Julich- 
Bonn group produce different results as compared to 
those performed with the corresponding parameters of 
the Nijmegen group. The parameters of model A of 
Ref. [p5| provides the best agreement with the data. 

An obvious extension of the present work is to calcu- 
late the cross sections for the pp — > pY}K + reaction, for 
which the measurements have recently been performed 
at COSY 0. This will also lead to the inclusive K+ 
cross sections in the elementary nucleon-nucleon colli- 
sions which are the necessary input to the transport 
model calculations of the strangeness production in the 
heavy ion collisions. This work is currently underway by 
extending the model to include the excitation of delta 
isobars A(1910) and A(1920) which are four star and 
three star resonances respectively. Since, the on-shell 
Y*(1650) — > YiK decay is not allowed, techniques similar 
to those described in section (D.ii) will have to be used 
to calculated the coupling constant for this vertex. It 
would also be interesting to calculate the invariant mass 
spectrum of the hyperon K + pair which is expected to 
provide further information about the various resonance 
contributions to this reaction. This will be reported in 
a future publication. Extensions of the present theory 
to incorporate the unitarity perhaps on the lines of the 
if-matrix approximation |57j is also desirable. 
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